The frustrated Brownian motion of nonlocal solitary waves 
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We investigate the evolution of solitary waves in a nonlocal medium in the presence of disorder. 
By using a perturbational approach, we show that an increasing degree of nonlocality may largely 
hamper the Brownian motion of self-trapped wave-packets. The result is vahd for any kind of 
nonlocality and in the presence of non-paraxial effects. Analytical predictions are compared with 
numerical simulations based on stochastic partial differential equations. 



Wave-packets may display particle-like behavior in the 
presence of non-linearity. Solitary waves (SW) and soli- 
tons are the non spreading solutions of the relevant non- 
linear wave equations that describe such a phenomenon. 
These self-trapped beams have been observed in a vari- 
ety of physical systems, ranging from oceanic waves to 
Bose Einstein condensates (BEC) [l|, l2|. Over the years 
the role of a nonlocal nonlinear response, with special 
emphasis on the optical spatial solitons (OSS) [3|, ap- 
peared with an increasing degree of importance |4l4lC| ; on 
one hand because it must be taken into account for the 
quantitative description of experiments and, on the other 
hand, because it is a leading mechanism for stabilizing 
multidimensional solitons ll|. Nonlocality in nonlinear 
wave propagation is found in those physical systems ex- 
hibiting long range correlations, like nematic liquid crys- 
tals (LC) [5|, photorefractive media (PR) [l^l, thermal 
[a, llSl, llJ] and thermo-difFusive tl5| nonlinear susceptibil- 
ities, soft-colloidal matter (SM) [1^, BEC [13, 111, and 
plasma-physics 19l |20| . When nonlocality is known to 
play a role, randomness is typically present and results 
from the unavoidable material fluctuations, and its role is 
fundamental for all-optical logic gates and soliton driven 
devices [JlliJl. 

In recent years, widespread investigations dealt with 
the interplay between randomness and nonlinearity, with 
emphasis on Anderson localization and SW |23l429l | . Un- 
derstanding the interplay between nonlocality and ran- 
domness is hence a fundamental subject in the theory of 
self-trapped waves. However, while the literature on soli- 
tons in random media is vast (see, e.g., [30| and references 
therein), stochastic nonlocal models were considered only 
in a few works [3l|, l32] and a quantitative theory for the 
SW fluctuations was not reported. 

In this Letter, we consider the effect of randomness 
on the stochastic dynamics of a nonlocal SW and show 
that, as the degree of nonlocality increases, the amount 
of fluctuations diminishes, and vanishes for an infinite- 
range nonlocal response. This is somehow a counter- 
intuitive result because in such media the nonlinear per- 
turbation extends far beyond the spatial extension of the 
self-trapped wave, which is hence expected to be influ- 
enced by the material fluctuations on a scale much larger 



than in the local case. 

The model we consider is written as 

Q{p)^\i^?+ri{x,t) 



(1) 
(2) 



where ip is the relevant wave fleld, x is the position 
vector, Q in a, linear differential operator, which does 
not include derivatives with respect to the evolution co- 
ordinate t, and ri{x,t) is a Langevin noise such that 
{ri{x,t)r]{x' ,t')) = rilj6{x-x')6(t—t'). The origin of ry de- 
pends on the speciflc physical problem: (i) temperature 
(nematic director) p fluctuations for thermal (LC) media; 
(ii) SM particle density p fluctuations; (iii) space-charge 
fleld p fluctuations for PR (eventually induced by mod- 
ulation of the background field) ; (iv) finite-temperature 
results in terms like 77 for plasmas and BEC [33| . 

In the Fourier domain ([2]) is written as p = >S'(g)(|^p -I- 
fj), where the tilde denotes the Fourier transform, S{q) is 
the "structure factor," [16[ , and the corresponding Green 
function is denoted by K(x), such that 



idt^p + V^V + Vix, t)tl) + i>K* \ip\^ = 0, 



(3) 



where V{x, t) — K * 1^ in a, colored random noise (the 
asterisk "*" denoting the x— convolution integral). The 
local regime corresponds to K{x) = 5{x)^ while in 
the highly nonlocal regime K{x) = Kq [J|. Let (j) — 
tpexp{il3t), Eq.([3l) is written as 

idt(l) + Vl<j)-l3(t) + (l)K*\(l)\^ = is(x,^p,i;a:,ip^a:,...,t) (4) 

where s is taken as a perturbation term, depending on 
-0 and its transverse derivative at any order, and /3 is 
the nonlinear wave- vector. Eq.(l4]) is a generalization of 
([3|), accounting for any kind of perturbation, e.g., in the 
presence of material losses s = —acj), with a the loss 
coefficient. Eq. (|3]) corresponds to s = Vcj). 

Soliton perturbation theory was previously developed 
for one dimensional (ID) solitons of the integrable local 



nonlinear Schroedinger (NLS) equation (see, e.g., |34|). 
and it is based on the knowledge of the exact soliton 
solutions. Here this approach is generalized to a non- 
integrable model, in the presence of an arbitrary non- 



locality, and then applied to the SW Brownian mo- 
tion. These are given, in the absence of the perturbation 
(s = 0), by the real valued solutions u{x) of 

- I3u + Vlu + u I K{x - x')u^{x')dx' = 0. (5) 



In absence of perturbation, the general solution of (j4]) is 
dependent on /? and on a given number of parameters, 
due to the Lie symmetries of the model. To simplify 
the notation, we consider hereafter the ID case, while 
the results below equally apply to the multi-dimensional 
case. The un-perturbed solution is written as 

(j)o^u{x- X + 2nt, P) exp(i6' - iVlx - iVt^t) (6) 

where translational, gauge and Galilean invariance are 
taken into account: X is the center of the self-trapped 
wave, 9 is the phase, and £7 is the momentum. Due to 
the Galilean invariance, the analysis can be limited to 
SW with zero velocity (fJ = 0). By letting (j) = (j)^ -\- cjji^ 
the linearized evolution equation is 



with 



dt4)i = C{(t)i) + s, 



C{(j)l) = -iPcf)! + i(j)l,xx + 

i(j)oK * {(t)Q(l3l + 0*</)i) -I- i(j)iK * |0o|^ 



(7) 



(8) 



By introducing the scalar product (a, 6) — ^ J a*bdx, 
it can be readily verified that C satisfies the following 
relation (a, C{b)) = (/3(a), b) with C the adjoint operator 



£(01) = z/3(/.i - 
K * (0o(/)t - 05( 



i(t>i,xx - i^ox 
)i) ~i(t)iK * |0oP 



(9) 



and such that C{ia) = —iC{a). Without loss of gen- 
erality, the first order perturbation can be decomposed 
in a term representing a small variation of the solitary- 
wave parameters and the remaining part, denoted as the 
radiation term (jjr. The former is proportional to the 
derivatives of 0o with respect to the various parameters, 
X ,Q,9,/3 , and the expression for (jji is written as 

(/-i = fx6X + fe&e + fp6/3 + {fn - Xfg)6n + 0„ (10) 

while having introduced the auxiliary functions 



fe = i4>o 

fx = dx<l>o 

fn = -i{x - X)(l)o, 



(11) 



and 5X{t), 5e{t), 6p{t) and 5Vl{t) being the time- 
dependent perturbations to the bound state parameters. 
Any of the auxiliary function / is such that dtf — £(/), 
the adjoint functions / are defined by 9t(/, /) = 0. They 
are given by fg = ifpjp = -ife, fa = -ifx, fx = ifn 



and are such that {fa,fb) = J^a^aM with a and b two 
symbols in the ensemble {X,n,9,f3). By direct integra- 
tion it turns out that Afe = Afp = {l/2){dP/df3) and 
Afx = A/o = (1/2)P, with P is the propagation in- 
variant SW norm, or power, P = J'|(/)op(ix = P{(3). 
dP/dp = P' > due to nonlocal SW stability Q. Fol- 
lowing the original argument in [35| , the functions / and 
/ are localized around the SW position X, as a conse- 
quence, their scalar products with 0^ vanish because the 
radiation term spreads for long times, and (/, (pr) = is 
invariant. This assumption holds true at the lowest or- 
der of approximation in s, as based on the integrable NLS 
equation, and is confirmed in the non-integrable case con- 
sidered here by the agreement with numerical simulations 
reported below. 

By using the previous expressions in ([7]) and project- 
ing on the adjoint functions, the following equations are 
derived for the dynamics of the SW parameters 



XSn = (5/3 



SX = -25rj 



2^x 
P 



2^ 
P' ' 



6(3 






25n 
P 



(12) 



where Sa = (/q,s), and the dot is the i— derivative. 

Eqs. (J12p hold for any s; for a random perturbation in 
the density p, as given by ri{x, t) above, they become 

59 = X5^ + '^^+^7^5/"'^'^" ^^^^ * •^^'^^' 

6/3 = 0, SX ^ -2n, (13) 

6n = I u{x- X)ux{x - X){K * f)dx, 

from which 

sn{t) = -^J^ 

J Ju{x- X)ux{x - X)K{x - x')f{x', t')dx'dxdt'. 

(14) 
By writing 5Vl{t)6Q.{t') after P^ and averaging over dis- 
order leads to 

{5^{t)5D.{t')) = ^C min(f , t') (15) 



where 

C = / / / u{xi - X)u{x2 - X)u{x3 -X)x 
Ux{x3 — X)K{xi — X2)K{xg, — X2)dxidx2dx3 

The deviation from the mean position is found as 

5X{t) = -2 [ 5n{t')dt' 
Jo 

from which {5X) ~ and 



(16) 



(17) 



{sx{ty') = 4 / / {n{h)n{t2))dhdh - ^?^t" 



3P2 



(18) 



The previous result is the nonlocal counterpart of the so- 
called Gordon-Haus effect, first introduced for describing 
the random fluctuations of solitons in amplified light- 
wave systems [3J, |35|- Ea. ([T51) states that the random 
fluctuations, measured by (5X(i)^), grow with the cubic 
of the propagation distance, and decay with the square 
of the SW power. A significant role is played by quantity 
C, which, in general, depends on the specific soliton and 
nonlocality profile. In the local regime iir(2:—y) — S{x~y) 
it is 



C 



[u{xi — X)ux{xi — X)] dxi. 



(19) 



The relevant result is found in the highly nonlocal limit 
K{x) — )► Kq, which gives, for a bell-shaped soliton profile 

[u{x) =u{-x)], 




C 



^0 



u{xi — X)ux{xi — X)dx\ 



0, 



(20) 



irrespective of the specific shape of K{x). As a result, in 
the highly nonlocal regime the random, fluctuations of the 
fundamental soliton vanish. Physically, this corresponds 
to the fact that spectral power density of the noise is av- 
eraged out by a narrow S{q) as the degree of nonlocality 
increases. We stress that this results is independent of 
the specific kind of nonlocality. For example, with refer- 
ence to the exponential nonlocality S{q) — (1 



aV) 



11| , we show in figure [T] the C parameter versus ct^ for 
various P (note that /3 changes along each curve, because 
cr^ varies) as calculated after the bound-state solutions 
of Kq.©. As expected when a^ increases (ct^ = cor- 
responds to the local case), the predicted fluctuation de- 
creases. Eq.([5]) is also valid in the two-dimensional case 
for each transverse coordinate. 

To validate the previous analytical results, we resorted 
to the numerical integration of the stochastic partial dif- 
ferential equation resulting from a ID exponential non- 
locality; we adopted a pseudo-spectral stochastic Runge- 
Kutta method 36|, |37| . Figure [2^ shows a typical evolu- 
tion starting from a bound state and displaying the ran- 
dom deviation of the SW. In flgure[2j3, we report various 
trajectories for a fixed SW power. Figure |3] shows the cal- 
culated standard standard deviation for various degrees 
of nonlocality: the analytical prediction after Eq.(fT8|) is 
indistinguishable from the numerical results. 

Before concluding, we consider the effect of non- 
par axiality on the SW fiuctuations. Ultra-thin nonlo- 
cal OSS were considered in 16|, |38[; in this framework. 



non-paraxiality at the lowest order is described by the 
perturbation s = —e'!/)^;^;^;^;, where e is ratio between the 
wavelength and the spatial beam waist (within some nu- 
merical constants, see 16|). Such a term is orthogonal to 
all the adjoint functions, with the exception of fg] this im- 
plies that Ea. ((T5)) for the solitary- wave fluctuations still 
holds true in the ultra-focused regime, with the addition 



FIG. 1: (Color online) Amplitude of the fluctuations C for the 
exponential nonlocality versus cr^, after Ea. (|18p . for various 
soliton powers P. 
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FIG. 2: (a) Typical dynamics of a solitary wave [Eq.®] as ob- 
tained from the numerical solutions of the stochastic equation 
(Q in the presence of randomness and exponential nonlocality 
(P = Q, a^ = 5, r)N ~ 0.01); (b) Center of mass trajectories 



of the same bound-state 
realizations (r]N ~ 0.01). 



(P = 6, a ^ 5) for 100 disorder 



of a linear increase of the soliton phase along propaga- 
tion: 



9 -^t^ 

^non-paraxial — n/ 



UpUx 



r.dx. 



(21) 



corresponding to the perturbation to the nonlinear wave- 
vector of the bound-state due to the non-paraxial term. 
Conclusions — We have theoretically shown that nonlo- 
cality largely affects the dynamics of a solitary-wave in 
the presence of disorder; this turns out into a random 
walk of the self-trapped beam position, which is ham- 
pered by the filtering action of the nonlocal response, 
and ideally vanishes for a infinite degree of nonlocality. 
These results are expected to be specifically relevant for 
plasma-physics, Bose-Einstein condensates, thermal and 
thermo-diffusive media, liquid crystals and soft-colloidal 
matter, and suggest to employ highly nonlocal media for 
routing information by solitons in order to moderate the 
effect of randomness. 
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FIG. 3: (Color online) Comparison of the numerically (con- 
tinuous lines, P = 6, tjn ~ 0.01) and theoretically (circles for 
a^ = 2, squares for a^ — 5 and diamonds for a^ = 10) calcu- 
lated standard deviation of the solitary wave position versus 
the evolution coordinate for three degrees of nonlocality and 
100 disorder realizations. 
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